Surface bound states in the continuum 
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We introduce a novel concept of surface bound states in the continuum, i.e. surface modes 
embedded into the linear spectral band of a discrete lattice. We suggest an efficient method for 
creating such surface modes and the local bounded potential necessary to support such embedded 
modes. We demonstrate that the embedded modes are structural stable, and the position of their 
eigenvalues inside the spectral band can be tuned continuously by adding weak nonlinearity. 
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Soon after the emergence of quantum mechanics, von 
Neumann and Wigner suggested [T] that certain poten- 
tials could support spatially localized states within the 
continuum spectrum, i.e. bound states with the ener- 
gies above the potential barriers. Since 1929, when this 
remarkable proposal was published, the bound states in 
the continuum were regarded as a mathematical curiosity, 
even though such potentials were suggested to occur in 
certain atomic and molecular systems [H |3] . The subse- 
quent experiments with semiconductor heterostructures 
provided the direct observation of electronic bound states 
above a potential well localized by Bragg reflections [I]. 

In addition to the physics of electronic structures and 
quantum dot systems, this topic attracted a lot of atten- 
tion in optics [SHI] i where it was very recently shown that 
the optical bound states can be generated in an optical 
waveguide array by decoupling from the continuum by 
virtue of symmetry only [8]. 

In this Letter, we extent the pioneering concept of 
von Neumann and Wigner pQ into two novel directions. 
First, we demonstrate that the bound states can exist in 
systems of a semi-infinite extent as surface bound states 
in the continuum. We suggest and implement a novel 
method for creating square-integrable, discrete surface 
modes embedded into a linear spectrum. Such surface 
modes can be regarded as a novel type of localized sur- 
face Tamm-like states with energies in the continuum (i.e. 
"embedded Tamm modes" ) . Second, we study the prop- 
erties of such embedded states in the presence of nonlin- 
earity and demonstrate that the mode location inside the 
band can be tuned continuously by changing the mode 
amplitude. Importantly, our modes appear in entirely 
asymmetric systems and therefore, cannot be reduced to 
the bound states in infinite systems where the decoupling 
from the continuum occurs due to the symmetry condi- 
tions. 

We consider a semi-infinite, one-dimensional linear lat- 
tice, in the presence of a site energy distribution {e„} 
(Fig[T]). In optics, this could correspond to a set of 
weakly-coupled optical waveguides [5] , each of them char- 
acterized by a propagation constant e n and centered at 
x n = na. In the coupled-mode approach, we expand 
the electric field E(x, z) as a superposition of the fun- 
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FIG. 1. (Color online) Example of the site energy distribu- 
tion (grey bars) and its associated surface bound state (solid 
curve) with energy inside the continuum spectrum. 



damental modes centered at each waveguide, E(x, z) = 
J2 n Cn(z)4>{x — na ), where <fi(x) is the waveguide mode. 
We pose C n {z) = C n exp(iAz). After inserting this into 
the paraxial wave equation, one obtains the stationary 
equations for the mode amplitudes 



(-A + e n )C n + V{C n+1 + C n _i) =0, n > 1, 
and the similar equation for the edge waveguide, 

(-A + ei)Ci + V C 2 = 0, n=l. 
From these equations, we can formally express 

A ]'{'--- ) if n = l; 



\ — V 




if n > 1. 



For a homogeneous system, we take e„ = 0, and 
C n = sin(fcn) and A = 2V r cos(fc). 
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For e„ 7^ 0, and in the spirit of the concept of von Neu- 
mann and Wigner [J, we look for a localized surface 
mode in the shape of a modulated wave of the form 



C n = sin(nfc) /„ 



(5) 
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where / n — ^ for n —> oo, in order to realize a localized 
state. After inserting this ansatz into ([lj, we obtain 

[ A - 2V(f 2 /fi) cos(fc) ifn = l 

e„ = < A - V"{ (/ n+ i// n )[cos(fc) + sin(fe) cot (fen)] 

[ +(/n-i//n)[cos(fe) - sin(fc) cot(fcn)] } if n> 1. 

(6) 

and we have Iim n _ J , o0 e„ = 0, provided that 
lim„^ 00 (/ n+ i//„) = 1. 
Let us take 



Jn 



+ 1 



J 71 



= (1 - S n ) 



(7) 



where 5 n < 1. From this, we can solve formally for /„: 



fn = J] - 5 ™) 



(8) 



which can be rewriten as 
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/„ = cxp < iog(i - s m ) 



(9) 



In the limit n — > oo, and using that 8 m < 1, we can 
approximate this by 



FIG. 2. (a) Discrete envelope function /„ vs. n. (b) Ratio of 
the envelope functions at the neighboring sites, (c) Embedded 
mode profile, (d) Close up of the panel (c). (e) Site energy 
distribution and (f) participation ratio of all eigenvectors of 
the linear modes, (k = 0.56, N = 533). 



/oo sa exp < - ^ 



(10) 



where, we want — > 0. This implies J2m=i $m = 00 • 
A good trial function for 5 n is 

S n = -4= sin 2 (nfc) sin 2 ((n + l)jfe). (11) 

The presence of the sine terms is not accidental; we need 
them to counteract the presence of the two cot(nfc) terms 
in ([6| that may otherwise lead to possible divergences. 
In this way, we get a smoother site energy distribution. 
A useful parameter to quantify the degree of localiza- 
tion of a state, is its participation ratio R, defined by, 
R = (En |C„| 2 )7£„ |C„| 4 . For localized modes, R « 1 
while for extended states R « JV, where N is the number 
of sites in the lattice. 

Figure [2] shows results for a lattice of N = 533 
sites, using the trial function (11) and k = 0.56. The 
mode approaches zero slowly but surely. The asymp- 
totic decay of the envelope at large n values can be es- 
timated, using the Euler-Maclaurin formula to be C n ~ 
sin(fen) exp(— a(k)y/n), with a(k) — (2+cos(2fc))/4. Fig- 
ure [2] also shows the site energy distribution e„ and par- 
ticipation ratio R of all modes inside the band (out- 
side the band, there are 10 "impurity" localized states), 
and we see that our candidate for embedded mode (with 
eigenvalue A = 1.695) has the lowest of them all, R ps 3. 
The next higher R value is w 85. Figure [3] shows the 
states inside the band that are closest in energy to the 



embedded state. The embedded state is the only state in- 
side the band whose amplitude decreases to zero at large 
distance from the surface (n — 1), while all the rest of 
the band states are extended. 

Next, we proceed to check the structural stability of 
the embedded mode. That is, whether the mode is sta- 
ble against perturbations caused, for instance by errors in 
the form of the site energy distribution (refraction index 
distribution in optics). This can happen during an at- 
tempt to reproduce experimentally the embedded mode. 
We took a lattice of 333 sites and examined two cases: 

(a) The energy site distribution {e„} is replaced with an- 
other e„ —¥ e„ + <5„, where d n is a random number taken 
from a uniform random distribution whose width is pro- 
portional to the e„ at a given site. For instance, we took 
S n £ [ -0.1 |e° W |,0.1 |e° w | ]. The old and new energy site 
distribution look nearly the same, and as a result, the old 
given A, just shifts a little bit. The number of states out- 
side the band does not change and we still have an embed- 
ded mode, surrounded by extended modes, as in Fig. 3. 

(b) Next, we try a more drastic perturbation, where 5 n 
is now drawn from a uniform random distribution which 
does not depend on site position: e™ 6 ™ = e° w + £„ where, 
5 G [—0.1, 0.1]. We see in this case that, even though the 
e n becomes significantly distorted far from the boundary, 
the embedded state, with eigenvalue 1.89208 survives, 
with a different eigenvalue 1.85174. Whether the new 
eigenvalue is smaller or large than the original one, de- 
pends on the random realization. Figure [4] shows the old 
and new site energy distribution, while Fig. [5] shows the 



3 





0.1 




0. 




-0.1 








0.1 



-0.1 





0.1 




0. - 




-0.1 




0.1 






0.1 




0. 






0. 




0.1 












-0.1 






-0.1 





lllll 01 




o. 




1 -o.i 




0.1 
0. 





0.1 








-0.1 





300 



3(11) 




3()() 



FIG. 3. States in the spectrum band that are the closest in 
energy to the embedded mode (middle panel) (N = 333, k = 
0.33). 



FIG. 5. Perturbed states in the continuum spectrum band 
that are the closest in energy to the embedded mode (left 
middle panel) (N = 333, k = 0.33). 



new band states surrounding the new embedded mode. 
The spatial profile of all of them maintain their extended 
nature, save for a small tendency towards localization, in 
agreement with Anderson localization theory. As long as 
the disorder is small and Anderson's localization length 
is much larger than the dimensions of the lattice, the 
localized embedded state is well defined. 

We can then conclude that the embedded state is struc- 
turally stable against small perturbations. 

Finally, we address the issue of possible tunability of 
the embedded state. Staring from a given embedded 
mode A, we would like to be able to change its position 
inside the band, without altering the original site energy 
distribution. An attractive way to do this, accessible in 
optics, is by introducing a small amount of nonlinearity 
into the system. For Kerr nonlinearity, all state eigenval- 
ues will be shifted by the amount 7 |C„| 2 . For extended 
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FIG. 4. Original (left) and randomly perturbed (right) site 
energy distribution (N = 333, k = 0.33). 



states, the effect will be negligible, and the only state 
that can be affected, is the localized one. In the presence 
of a nonlinear response, the discrete equations Eqs. ([!]), 
pi) become: 



(12) 



(-A+en+^CnnCn+VtCn+i+CVi) = 0, 71 > 1, 

and the equation for the edge waveguide 

(-A + e 1 + 7 |C 1 | 2 )C 1 + V C 2 = 0, n=l, (13) 

where 7 = 1(— 1) denotes attractive (or repulsive) non- 
linearity. At this point, it is useful to make the change 
of variables (f>„ — C n j\/P, where P = J2 n |C«.| 2 is 
the total power. The effective nonlinearity parameter 
is now x = 7-P, and the <f) n are normalized to unity: 
l^nl 2 = 1- The idea is to start from the linear em- 
bedded state with given eigenvector A at x = 0. Then, we 
gradually increase or decrease \ an d follow the evolution 
of its eigenvalue and spatial profile, by solving Eqs.(fT2l, 



(131 in a self-consistent manner. 

Results from this procedure are shown in Fig. [6] We 
see that the eigenvalue of the embedded state can indeed 
be tuned to occur at any value inside the band, by means 
of a small amount of focussing or defocusing nonlinearity. 
The embedded state profile does not change perceptible 
during this process, as expected from the above discus- 
sion. 

The procedure described above for one-dimensional 
semi-infinite discrete lattices could be generalized to 
higher dimensions. The simplest of such cases is a semi- 
infinite square lattice where the surface modes can ex- 
ist near the edges and corners of a large square lattice. 
There, and given the separability of the Hamiltonian that 
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FIG. 6. Eigenvalue of the embedded surface mode as a func- 
tion of the nonlinearity strength, for N = 133 and k = 0.33. 
Black circle marks the position of the embedded eigenvalue in 
the absence of nonlinearity. 

gives rise to the evolution equations, the spatial profile 
of the embedded mode can be written as C n>m — 4>n4>m, 
where (f> n = sm(nk x )f% and <p m — sm(mfcy) where 



the envelope functions and f% are the same. The 
eigenvalue of the embedded mode is A = X x + X y and the 
energy site distribution is simply e n<m = e n + e m , each is 
given by the appropriate modification of Eq. Q . 

In summary, we have suggested and studied a novel 
type of bound states localized at the edge of a semi- 
infinite discrete lattice with the eigenvalues embedded 
in the continuous spectrum. We have demonstrated a 
procedure to generate square-integrable, surface local- 
ized modes embedded in the continuum, as well as the 
site energy distributions needed to produce such modes. 
We have shown that these new embedded modes are 
structurally stable, and their location inside the band 
can be tuned by weak nonlinearity. We believe the idea 
demonstrated here may be useful in other fields, includ- 
ing atomic systems, quantum-confined structures, as well 
as a variety of photonic structures. 

The authors acknowledge support from FONDECYT 
Grant 1080374 and Programa de Financiamiento Basal 
de CONICYT (FB0824/2008), and the Australian Re- 
search Council. 



[1] J. von Neumann and E. Wigner, Phys. Z. 30, 465 (1929). 
[2] F. H. Stillinger, J. Chem. Phys. 45, 3623 (1966); F. H. 

Stillinger and T. A. Weber, Phys. Rev. A 10, 1122 (1974); 

F. H. Stillinger and D. E. Stillinger, Phys. Rev. A 10, 1109 

(1974) ; F. H. Stillinger and D. R. Herrick, Phys. Rev. A 
11, 446 (1975). 

[3] A. K. Jain and C. S. Shastry, Phys. Rev. A 12, 2237 

(1975) . 

[4] F. Capasso, C. Sirtori, J. Faist, D. L. Sivco, S-N. G. Chu, 
and A. Y. Cho, Nature 358, 565 (1992). 



[5] D.C. Marinica, A.G. Borisov, and S.V. Shabanov, Phys. 

Rev. Lett. 100, 183902 (2008). 
[6] E.N. Bulgakov and A.F. Sadreev, Phys. Rev. B 78, 075105 

(2008). 

[7] N. Moiseyev, Phys. Rev. Lett. 102, 167404 (2009). 

[8] Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, 

A. Szameit, and M. Segev, Phys. Rev. Lett. 107, 183901 

(2011). 



